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ABSTRACT 



We used the Very Long Baseline Array (VLBA) and the European VLBI 
Network (EVN) to perform phase-referenced VLBI observations of the three 
most powerful maser transitions associated with the high-mass star-forming re- 
gion G28.87+0.07: the 22.2 GHz H 2 0, 6.7 GHz CH 3 OH, and 1.665 GHz OH 
lines. We also performed VLA observations of the radio continuum emission at 
1.3 and 3.6 cm and Subaru observations of the continuum emission at 24.5 /zm. 
Two centimeter continuum sources are detected and one of them (named "HMC" ) 
is compact and placed at the center of the observed distribution of H 2 0, CH 3 OH 
and OH masers. The bipolar distribution of line-of-sight (l.o.s) velocities and 
the pattern of the proper motions suggest that the water masers are driven by 
a (proto) stellar jet interacting with the dense circumstellar gas. The same jet 
could both excite the centimeter continuum source named "HMC" (interpreted as 
free-free emission from shocked gas) and power the molecular outflow observed 
at larger scales - although one cannot exclude that the free-free continuum is 
rather originating from a hypercompact H II region. At 24.5 /im, we identify 
two objects separated along the north-south direction, whose absolute positions 
agree with those of the two VLA continuum sources. We establish that ~90% 
of the luminosity of the region (~2xl0 5 L Q ) is coming from the radio source 
"HMC" , which confirms the existence of an embedded massive young stellar ob- 
ject (MYSO) exciting the masers and possibly still undergoing heavy accretion 
from the surrounding envelope. 

Subject headings: ISM: individual objects (G28.87+0.07) - ISM: kinematics and 
dynamics - masers - techniques: interferometric 
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INTRODUCTION 



Understanding the process of high-mass star formation represents a challenge from both 
a theoretical and observational point of view. By now, three main com peting concepts of 
massive star formation have been discussed in the recent literature (see 



Zinnecker fe Yorke 



20071 . and references therein): f ormation throu 



gh coalescence of lower-mass stars in 



extremely dense stellar clust ers ( 



Bonnell et al. 



turbulence supported cores ( jMckee fc Tan 



19981 ); monolithic collapse of isolated, 



20031 ); and competitive accretion where the 



most massive protostars benefit of the gravitational well c reated by the lower-m ass stars 



Bonnell fc Bate 



20061 ). The 



surrounding them, thus boosting the accretion rate (e.g. 
available observational evidence is still inadequate to discern between these different 
scenarios, mostly because observations of massive young stellar objects (MYSOs) are 
complicated by the large distances of typically several kpc, implying small angular scales, 
and the fact that such objects are born still deeply enshrouded in dense dusty envelopes, 
opaque to optical and near-infrared wavelength radiation. 

A possible way to distinguish between different models is to explore the physical 
conditions and the dynamical processes of the gas close to the protostars. However, present 
(sub)millimeter interferometers and infrared (IR) telescopes attain angular resolutions 
>0'/25, corresponding to thousands of astronomical units (AU) at the typical distances of 
MYSOs of several kpc. This is insufficient to investigate the kinematics of the gas in the 
immediate surroundings of the forming star, where crucial phenomena such as accretion, 
rotation, and ejection are likely most prominent. The above mentioned limitation can be 
bypassed by Very Long Baseline Interferometry (VLBI) observations of maser lines, which 
attain linear resolutions of a few AU. 

A few years ago, we started an observational project to study the high-mass 
star-forming process by comparing interferometric images of thermal lines of molecular 
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tracers with multi-epoch VLBI studies in three well known maser species (OH, CH3OH 
and H 2 0). At present, three sources (IRAS 20126+4104, G16.59-0.05 a nd G23.01-0.41) 



have been already analyzed (jMoscadelli et al 



2011 



Sanna et al 



2010a 



bj). The results 



demonstrate that the synergy between VLBI, multi-species maser observations and 
(sub)mm interferometric observations of thermal molecular lines is crucial to achieve a 
multi-scale picture of the environment of newly formed massive (proto) stars. In the present 
article we report on our VLBI maser study of the high-mass star-forming region (HMSFR) 
G28.87+0.07, and compare our results with available complementary interferometric data 
of continuum and thermally excited molecular line emission, as well as with sub-arcsec 
resolution mid-IR images obtained with the Cooled Mid-Infrared Camera and Spectrometer 
(COMICS) on the Subaru telescope. 

In Section 2, we provide a review of previous observations toward this region. Section 3 
describes the VLBI observations of the maser transitions, the Very Large Array (VLA) 
observations of the radio continuum emission at 1.3 cm and 3.6 cm, and the mid-IR 
observations performed with the Subaru/COMICS telescope. In Section 4, we present 
the distributions of the three types of masers, the water maser proper motions, and the 
properties of the radio continuum. The results are discussed in Section 5 and the conclusions 
are drawn in Section 6. 



2. THE HMSFR G28.87+0.07 



The H MSFR G28.87- 
of 7.4 kpcE ( Codella et all 



p. 07 ( alias IRAS 18 



1997 



Furuya et al. 



11—0338) is located at a kinematic distance 



20081 ) and is associated with a pc-scale 



1 Since the a mbiguity between th e near (ps6 kpc) and the far (^9 kpc) kinematical distance 



is not resolved (jSewilo et al. 



20041 ). we adopt the distance to the tangent point. 
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clump with a mass in exce s s of 10 3 Mp, and a b olometric luminosity of ~2xl0 5 L & 



(ILopez-Sepulcre et al. 



2009 



Faundez et al 



2004 ). From the observation of NH 3 , the 



systemic velocity of the re gion (K ys ) with res pect to the local standard of rest (LSR) is 



found to be 101.0 km s 1 ( Codella et al. 



19971 ). 



On a large scale, mid-IR observations at 8, 10.5, and 18.1 /jm have reso l ved the IRAS 



sourc e in two objects along the north-south direction (IDe Buizer et al 



2008 



2005 



Furuya et al. 



Figure 1). Bot h sources could b e associated with an Midcourse Space Experiment 



(MSX) point source (jEgan et al.l 120031 ) and a 2.2 /xm source with strong nea r-infrared 



exces s, lying 3"-4" to the southeast of the northern, brighter m id-IR object (ITesti et al 



19941 ) . Bolometer observations at 1.2 mm ( jFaundez et al. 



20041 ) have detected a peak 



of dust continuum emission associated with the two mid-IR sources and the near-IR 
source. Single-dish maps in the CS(5-4) and C 18 0(2-l) lines have revealed an arcmin-scale 



molecular clump and a monopolar (red-s 



(jShirley et al. 



2003 



Lopez- Sepulcre et al. 



lifted ) outflow is seen in the CO (2-1) line 



20091 ). Single- dish surveys 



rave detected emission 



in dense gas tracers as CS (7-6), CS(2-1), and CS(l-O) (IPlume et al. 



1992 



1996 



Anglada et al. 



Araya et al. 



19961 ). and in the HI 10a radio recombination line ( jSewilo et al. 



Bronfman 



et al. 



2004! : 



20071 ). which suggests the presence of a diffuse H II region. In addition, OH 
maser emission has been found in the main (1.665 and 1.667 GHz) hyperfme structure lines 
and in the 1.7 20 GHz satellite line, whereas the 1.612 GHz OH satellite has been observed 



in absorption (ILittle fc Cesarsky 



1982 



Caswell fc Havnes 



maser emission hay e been reported at 6.7 GHz (ICaswell et al. 



1983). W eak C H 3 OH and strong 



Genzel fc Downeslll979l ). respectively. 



19951 ) and 22.2 GHz 



Interfer ometric observation s at millimeter wavelengths have revealed multiple sources in 



this region (iFuruya et al. 



2008 



Figure 1) and show that the dust continuum peak coincides 



with the position of the northern, stronger mid-IR source. The mass (100 M ) from the 
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3 mm dust continuum emission and the rotational temperature (93 K) inferred from the 
CH 3 CN lines suggest that one is dealing with a hot molecular core (HMC). In this study, 
we will focus on this HMC, as it contains all the tracers of high-mass star formation that 
we have investigated with our obs ervations, namely th e ma sers and the fr e e-free and mid-IR 



continuum sources. As shown by 



Codella et al 



(119971 ) and 



Furuya et al 



(200 81), the HMC 



is seen in NH 3 , HCO + , and CH 3 CN line emission. The observations of ICodella et al. 



19971 ) 



were not sensitive enough to detect continuum emission at 1.3 cm, whereas our new images 
demonstrate the existence of a free- free continuum source inside the HMC (see below). A 
wide-angle, massive (>100 M ) bipolar outflow has been detected in the 12 CO(1-0) line, 



orien ted perpendicular to a CH 3 CN velocity gradient detected in the HMC (IFuruya et al. 



200 8]) . The absolu t e pos itions of the 1.665 GHz OH and 22.2 GHz H2O masers, derived by 



Forster fc Caswell! ( 119991 ) using the VLA with an accuracy of C/5, indicate that these maser 



emissions are associated with the HMC. 

All these features are consistent with the presence of a deeply embedded MYSO(s) 
inside the HMC, possibly accompanied by lower mass objects lying in the neighborhoods of 
the HMC itself. 



3. OBSERVATION AND DATA REDUCTION 

3.1. VLA: 1.3 cm and 3.6 cm Continuum 

We observed the HMSFR G28.87+0.07 (tracking center: R.A.(J2000) = 18 h 43 m 46.24 s 
and Dec.(J2000) = -03°35'30.4") with the VLA at X (3.6 cm) and K (1.3 cm) bands in 
both the C-array (in 2008 April) and A-array configurations (in 2008 November). At 3.6 cm 
in both configurations and at 1.3 cm in the C-array configuration, the continuum mode 
of the correlator was used, resulting in an effective bandwidth of 172 MHz. At 1.3 cm in 
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the A-array we used mode "4" of the correlator, with a pair of intermediate frequency (IF) 
bands, each with a width of 3.125 MHz (64 channels) centered on the strongest H 2 maser 
line and a pair of 25 MHz wide IF bands (8 channels each) sufficiently offset from the maser 
lines to obtain a measurement of the continuum emission. The two bands centered at the 
same frequency (detecting both circular polarizations) were averaged to decrease the noise. 

At X-band, 3C286 (5.2 Jy) and 3C48 (3.1 Jy) were used as primary flux calibrators, 
while 1832—105 (1.4 Jy) was the phase calibrator. For the K-band observations, the 
primary flux calibrator was 3C286 (2.5 Jy), the phase calibrator 1832—105 (1.0 Jy), and 
the bandpass calibrator (for the A-array only) 1733—130 (3.7 Jy). 

The data were calibrated with the NRAcJ^I Astronomical Image Processing System 
(AIPS) software package using standard procedures. Only for the A -array data at 1.3 cm, 
several cycles of self-calibration were applied to the strongest maser channel, and the 
resulting phase and amplitude corrections were e ventually transferred t o all the other 



line channels and to the K-band continuum data ( iMenten fc Reid 



19951 ) . This procedure 



resulted in a significant (at least a factor 2) improvement in the signal-to-noise ratio (SNR). 

Natural-weighted maps were made with task IMAGR of AIPS for both the continuum 
and the line data. Simultaneous observation of the line and continuum emission at K-band 
made it possible to obtain the relative position of the continuum image with respect to 
the H2O maser spots with great precision (~0'.'l). Finally, we identified the maser spots 
that appear in both our Very Long Baseline Array (VLBA) and VLA images and, using 
the accurate absolute position information obtained from the VLBA observations, we 
re-centered the VLA spots and continuum image. The correction thus applied turns out to 



2 The National Radio Astronomy Observatory is operated by Associated Universities, Inc., 
under cooperative agreement with the National Science Foundation. 
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be <20 mas in each coordinate. We conclude that the absolute astrometrical precision of 
the continuum maps presented in this paper is ~0'.'l. 



3.2. VLBI: Maser Lines 

We conducted VLBI observations of the H2O and CH3OH masers (at four epochs) 
and of the OH maser (at a single epoch) toward G28.87+0.07 at 22.2 GHz, 6.7 GHz, 
and 1.665 GHz, respectively. We give more technical details on the observations in the 
following subsections. In order to determine the masers' absolute positions, we performed 
phase-referencing observations in fast switching mode between the maser source and the 
calibrator J1834— 0301. This calibrator has an angular offset from the maser source of 2?4 
and belongs to the list of sources of the VLB A Calibrator Survey 3 (VCS3). Its absolute 
position is known with a precision better than ±1.5 mas, and its flux measured with the 
VLBA at S and X bands is 154 and 156 mJy beam -1 , respectively. Five fringe finders 
(J1642+3948; J1751+0939; J1800+3848; J2101+0341; J2253+1608) were observed for 
bandpass, single-band delay, and instrumental phase-offset calibration. 

The data were reduced with AIPS following the VLBI spectral line procedures. For 
a description of the general data calibration and the criteria used to identify individual 
masing clouds, derive their parameters (position, intensity, flux and size), and measure their 
(relative and absolute) prop er motions, we refer to the recent paper on VLBI observations of 



H 2 and CH 3 OH masers by 



Sanna et al. 



(j2010al ). For each VLBI epoch, H 2 maser absolute 



positions have been derived by employing two methods, direct phase-referencing and inverse 
phase-referencing, and the two procedures always gave consistent results. CH3OH and OH 
maser absolute positions have been derived only by using direct phase-referencing. The 
derived absolute proper motions have been corrected for the apparent proper motion due to 
the Earth's orbit around the Sun (parallax), the solar motion and the differential Galactic 
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rotation between our LSR and that of the maser source. We have adopted a flat Galactic 



rotation curve (Rq = 8.3 ± 0.23 kpc, O = 239 ± 7 km s 1 ) (IBrunthaler et al. 



Brunthaler et al. 


2011 


), and the 


5 +0 f 6 km s" 1 by 


Schonrich et al. 



(120101 ). who recently revised the Hipparcos satellite results. 



3.2.1. VLB A : 22. 2 GHz H 2 Masers 

We observed the HMSFR G28.87+0.07 (tracking center: R.A.(J2000) = 18 h 43 m 46.22 s 
and Dec.(J2000) = -03°35'29.9") with the VLBA in the 6 16 - 5 23 H 2 transition (rest 
frequency 22.235079 GHz). The observations (program code: BM244) consisted of 4 epochs: 
2006 April 23, 2006 June 30, 2006 September 28, and 2007 January 18. During a run of 
about 6 h per epoch, we recorded the dual circular polarization with 16 MHz bandwidth 
centered at an LSR velocity (Vlsr) of 101.0 km s _1 . The data were processed with the 
VLBA FX correlator in Socorro (New Mexico) using an averaging time of 1 s and 1024 
spectral channels. 

Images were produced with natural weighting, cleaned and restored with an elliptical 
Gaussian with full width at half maximum (FWHM) size of about 1.3 mas x 0.6 mas at 
a position angle (PA) of —4° (east of north), with small variation from epoch to epoch 
(shown in Table [2]). The interferometer instantaneous field of view was limited to ~2''7. In 
each observing epoch the on-source integration time was ~2.5 h resulting in an effective 
rms noise level of the channel maps (cr) in the range 0.006-0.015 Jy beam -1 . The spectral 
resolution was 0.2 km s -1 . 
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3. 2. 2. E VN: 6. 7 GHz CH 3 OH Masers 

We observed G28.87+0.07 (tracking center: R.A.(J2000) = 18 h 43 m 46.22 s and 
Dec.(J2000) = -03°35'29.9") with the European VLBI Network (EVN) in the 5 1 - 6 A + 
CH 3 OH transition (rest frequency 6.668519 GHz). Data were taken in 4 epochs (program 
code: EM061 and EM069) separated by ~1 year: 2006 February 28, 2007 March 18, 
2008 March 18, and 2009 March 17. In the first two epochs, the antennae involved in 
the observations were Cambridge, Jodrell2, Effelsberg, Hartebeesthoek, Medicina, Noto, 
Torun, and Westerbork. Since we found out that the longest baselines involving the 
Hartebeesthoek antenna (e.g., the Ef-Hh baseline is ~8042 km) heavily resolve the maser 
emission and do not produce fringe-fit solutions, the Hartebeesthoek antenna was replaced 
with the Onsala antenna in the third epoch. During a run of ~6 h, we recorded the dual 
circular polarization in two IF bands of width of 2 MHz and 16 MHz, both centered at 
an LSR velocity of 101.0 km s -1 . The 16 MHz band was used to increase the SNR of the 
weak continuum calibrator. The data were processed with the MKIV correlator at the Joint 
Institute for VLBI in Europe (JIVE-Dwingeloo, the Netherlands) using an averaging time 
of 1 s and 1024 spectral channels for each observing bandwidth. 

Images were produced with natural weighting, cleaned and restored with an elliptical 
Gaussian with FWHM size of about 12 mas x 6 mas at a PA of 21°, with small variations 
from epoch to epoch (shown in Table [2]). The interferometer instantaneous field of view 
was limited to ~9'/2. Using an on-source integration time of ~2.2 h, the effective rms noise 
level of the channel map (a) varied in the range 0.005-0.01 Jy beam -1 . The 2 MHz band 
spectral resolution was 0.09 km s _1 . 



3.2.3. VLB A: 1.665 GHz OH Masers 



We observed G28.87+0.07 (tracking center: R.A.(J2000) = 18 h 43 m 46.34 s and 
Dec.(J2000) = -03°35'29.9") with the VLB A in the 2 n 3/2 J = 3/2 OH transition (rest 
frequency 1.665401 GHz) on 2007 April 28 (program code: BM2440). During a run of ~6 h, 
we recorded the dual circular polarization with two IF bands of width of 1 MHz and 4 MHz, 
both centered at Vlsr = 103.0 km s _1 . The 4 MHz bandwidth was used to increase the 
SNR of the weak L-band signal of the continuum calibrator. The data were processed with 
the VLBA FX correlator in two correlation passes using 1024 and 128 spectral channels for 
the 1 MHz and 4 MHz bands, respectively. In each correlator pass, the data averaging time 
was 2 s. 

Images were produced with natural weighting, cleaned and restored with an elliptical 
Gaussian with FWHM size of 17.8 mas x 10.9 mas at a PA of 1°. The interferometer 
instantaneous field of view was limited to ~18'/5. The on-source integration time was 
~1.9 h, resulting in an effective rms noise level in each velocity channel of ~0.02 Jy beam -1 . 
The spectral resolution was 0.2 km s _1 . 



3.3. Subaru/COMICS: 24.5 |im Mid-infrared Continuum 



Kataza et al. 



20001) at the 



Using the mid-infrared imaging spectrometer (COMICS; 
Cassegrain focus of the 8.2 m Subaru Telescope, we carried out imaging observations of the 
24.5 /im emission toward G28.87+0.07 on 2008 July 15. For this purpose, we configured 
the Q24.5 filter, and employed chopping mode for subtracting the sky-background emission. 
The camera provides a field of view of ~ 42" x 32" with a pixel size o f O'/IS. Flux calibration 



was performed towards four standard sources listed in 



Cohen et al 



(119991): HD146051, 



HD186791, HD198542, and HD198542. We estimated the overall uncertainty in the flux 
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calibrations to be less than 10%. Astrometric calibration was performed by smoothing the 



24.5 (Ltm image to the same angular resq 



Galactic plane using MIPS: 



Carey et al 



ution as the 24 /j,m MIPSGAL (survey of the inner 



2009) image, i.e. 6", and overlaying the two images 



until a good match was found. Note that the heavy saturation of the MIPSGAL image 
does not hinder the comparison between it and the smoothed COMICS image. In practice, 
the emission in MIPSGAL is basically pointlike and the MIPS PSF presents a well defined 
circular structure even far from the (saturated) peak. The match between the smoothed 
COMICS image and the MIPSGAL image can thus be performed successfully on such 
circular features. Inspection by eye indicates a positional error of less than 1", comparable 
to the astrometric accuracy of MIPSGAL. 



4. RESULTS AND ANALYSIS 

This section reports our main results from the observation of the radio continuum, 
individual maser species and infrared emission towards the HMSFR G28.87+0.07. In the 
following, the term "spot" is used to refer to maser emission in a single velocity channel, 
whereas the term "feature" indicates a collection of spots emitting at a similar position 
over contiguous channels (i.e. an individual masing-cloud) . In this work, the parameters of 
maser emission (position, intensity, proper motion, etc..) are derived for maser features. 



4.1. 1.3 cm and 3.6 cm Continuum 

Our VLA measurements of the 1.3 and 3.6 cm continuum emission towards G28.87+0.07 
(in both the A- and C-array configurations) are summarized in Table [U, and shown in 
Figure [U At the angular resolution of the VLA C-array, the emission consists of two 
components, both detected at 1.3 and 3.6 cm. The northern and southern components, 
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separated by ~3'.'l, are labeled, respectively, "HMC" and "A" in Tabled] Th e former name 
i s due to the source being associated with the hot molecular core imaged by 



Furuya et al. 



( 120081 ). Source "A" is partly resolved with the C-array at 1.3 cm, and is completely resolved 



out with the A-array at both 1.3 and 3.6 cm. 

In the following, we focus our attention on the centimeter continuum source "HMC" , 
because this appears to be spatially associated with the observed 22 GHz H2O, 6.7 GHz 
CH3OH, and 1.665 GHz OH masers in G28. 87+0.07 (see Figures □ and H . While at 1.3 cm 
the most extended VLA A-array configuration partly resolves the emission of "HMC" 
(missing about half of the flux measured with the C-array), at 3.6 cm the emission is 
essentially compact in both configurations (yielding comparable integral fluxes). Note 
however that the VLA A-array image at 3.6 cm shows a spur protruding to the south. 
The peak positions of the 1.3 and 3.6 cm emissions coincide within the uncertainties. The 
absolute position of the A-array 1.3 cm image has been obtained by aligning positions of 
persistent water maser spots between the VLA and VLBA observations, and should be 
accurate to ~ O'.'l. 

Figure [2] shows the spectral energy distribution (SED) at centimeter wavelengths 
of "HMC". A spectral index (a) of 0.6 is derived from a linear fit of the VLA C-array 
1.3 and 3.6 cm and B-array 6 cm flux densities. The latter flux is an upper limit obtained 
from the GPS 6 cm Epoch 2 images downloadable from the MAGPIS web sitqj, which 
result to be the most se nsitive among the o ther GPS maps and the relevant map from the 



CORNISH survey I 4 ] (see 



a 



Purcell et al 



120081 ) . We also fit the SED at centimeter wavelengths 
with a model of a homogeneous H II region, with a Lyman continuum rate of 7.9 x 10 45 s~ x 
and a radius of 0'.'04. 



3 http:/ /third. ucllnl.org/gps 



4 available on http://www.ast.leeds.ac.uk/Cornish/public 
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We caution that the shape of the SED could be significantly affected by the different 
sensitivity to extended structures at different wavelengths. For the same VLA configuration, 
the shorter the wavelength, the larger can be the fraction of emission resolved out by the 
interferometer. Therefore, one should consider the possibility that the slope between 3.6 cm 
and 1.3 cm is steeper than 0.6. 



4.2. Masers Results 

4.2.1. 22.2 GHz Water Maser 



Forster fc Caswelll (119991 ) observed the HMSFR G28.87+0.07 with the VLA C-array 
(HPBW of 2'/2 x l'/6) and identified 13 water maser spots, mainly grouped in two clusters 
separated by ~0'.'2 in the northeast-southwest direction. The absolute position of the 
brightest spot (with a flux density of 61.3 Jy at the Vlsr of 106.3 km s _1 ) differs by ~0'.'l 
from that of the corresponding spot in our VLBA observation. The field of view of our 
VL BA observations c o vers t he whole area within which water maser emission was detected 
bv Forster fc Caswelll Jl999n . 



Figure [3] (upper panel) shows the total-power spectrum of the 22.2 GHz masers toward 
G28.87+0.07. This profile was obtained by averaging the total-power spectra of all VLBA 
antennae, after weighting each spectrum with the antenna system temperature (T sys ). 
Water maser absolute positions and LSR velocities are plotted in Figure HI Imaging the 
range of LSR velocities from 93.2 to 114.3 km s _1 , we detected 168 distinct water maser 
features distributed over a region of about 0'.'3 x 0'.'5. Most (92%) of them concentrate 
in a small area of about 0'.'3 x 0'.'2 centered on the VLA 1.3 cm continuum peak. The 
line-of-sight (l.o.s.) velocity distribution is remarkably bipolar, with red-shifted features 
located to the west of the radio continuum and blue-shifted ones to the east. While 
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the blue-shifted features mainly concentrate in two small (size <20 mas) clusters, the 
distribution of red-shifted masers consists of one narrow strip elongated east-west (size of 
falOO mas), together with two diffuse "clouds" of features spread to the north and to the 
south of the densely populated, narrow strip. It is also interesting to note the presence 
of two streamlines of water features (enclosed in dashed boxes in Figure H] upper panel), 
oriented north-south and extending over a few tenths of an arcsecond, on either side of the 
VLA 1.3 cm continuum. We took particular care in filtering out spurious emission owing 
to poor cleaning of the beam side lobes (elongated mostly along the north-south direction) 
and can confirm that these water maser streamlines are real. 

The individual feature properties are presented in Table [3j The position offsets are 
relative to feature #1. The absolute positions of feature #1 at the four observing epochs 
are given in Table HJ Feature intensities range from 0.1 to 54 Jy. The spread in LSR 
velocities ranges from 111.99 km s _1 for the most red-shifted feature (#107) to 93.33 km 
s _1 for the most blue-shifted one (#155). 

Only 61 features (36% of the total) persisted over at least 3 epochs, 34 of which lasting 
4 epochs. We calculated the geometric center (hereafter "center of motion" , identified with 
label #0 in Table |3j and indicated by a star in Figure |5]) of features with a stable spatial 
and spectral structure, persisting over the 4 observing epochs, and refer our measurement 
of proper motions to this point. With a time baseline of 9 months, relative (sky-projected) 
velocities of water maser features are derived with a mean uncertainty of ~23%. The 
magnitude of relative proper motions ranges from 1.9±2.8 km s -1 to 58.3±2.6 km s -1 , 
with a mean value of 24.0 km s -1 . Figure [5] shows the derived water maser proper motions. 
While the two blue-shifted clusters move to the east and northeast, most velocities of the 
red-shifted water masers point either to the south or to the west. In particular, looking at 
the distribution of the maser velocity vectors along the (east-west) narrow strip outlined by 
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the red-shifted features, we note that maser features closer to the "center of motion" move 
mainly to the south. Getting closer to the western end of the maser strip, the orientation 
of the maser velocity vectors changes gradually from south to west. The general pattern of 
water maser relative velocities appears to indicate expansion from a point on the sky close 
to the "center of motion" . Note also that at larger distances from the "center of motion" , 
maser velocities seem to be larger and better collimated along the (north)east-(south)west 
direction. In particular, we identify two H2O maser clusters located at the northeast and 
southwest ends of the maser distribution (named "Jf," and "J r ", respectively; see Figure |5]), 
which move with high velocities directed close to the northeast-southwest direction. 

Figure [6] (upper panel) reports the intrinsic absolute proper motions of the water maser 
features, derived from the observed absolute proper motion of feature #1 corrected for the 
apparent motion. The latter has been calc ulated adding the cont ribution of the parallax, 



the solar motion with respect to the LSR ( ISchonrich et al 



LO) 



2010), and the differential 



Galactic rotation, from a flat rotation curve with Galactic constants (Rp = 8.3 ± 0.23 kpc 



Brunthaler et al. 



( 1201 ll ). The obtained 



Go = 239 ± 7 km s _1 ) as recently determined by 
absolute proper motions differ from the relative proper motions presented in Figure |5] by a 
vector of amplitude ~40 km s -1 pointing to the southwest, which seems to be the dominant 
component of all proper motions. We think that this difference could be caused by the 
large systematic error due to the uncertainties in the source distance and solar and Galactic 
motion, and/or by a peculiar velocity of the maser source with respect to its LSR reference 
system. Our VLBA water maser observations were not designed to measure the source 
parallax. With a positional accuracy of ~0. 1-0.2 milli-arcseconds, and four epochs spanning 
(only) 9 months, the observed absolute proper motion of feature #1, not corrected for the 
apparent motion, looks similar to a straight line. From the small deviations from linear 
motion, we derive an upper limit of 0.2 mas for the amplitude of the sinusoidal parallax 
signature. This corresponds to a lower limit of 5 kpc for the source distance, in agreement 
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with the adopted kinematical distance of 7.4 kpc. 



We have checked the reliability of the relative, water maser proper motions by 
following an alternative approach. As described in Section. 14.2.21 we have also derived 



the absolute proper motion of the strongest 6.7 G 



VLBI observations of 6.7 GHz ( Sanna et al. 



2010a 



lz methanol ma ser feature. To date, 



2010 



b) and 12 GHz (IMoscadelli et al. 



2003 



20111 ) methanol masers toward several sources indicate that these masers usually 
show small proper motions, with typical velocities of only a few km s -1 . Assuming that 
also in G28.87+0.07 methanol masers move significantly slower than the water masers, we 
can use the absolute proper motion of the methanol masers to correct the absolute proper 
motion of the water masers, and derive the water maser velocities as approximately seen 
in a reference system comoving with the star. Figure [6] (lower panel) shows the absolute 
water maser proper motions after subtracting the absolute methanol maser proper motion 
(evaluated in Section. 14.2.2)) . Note that, within the measurement errors, the new absolute 
proper motions are fully consistent with the relative proper motions calculated with respect 
to the "center of motion" (Figure EJ). This result reinforces our assertion that the "center 
of motion" can define a suitable reference system to calculate water maser velocities. In the 
following discussion, we use the relative velocities of the water masers to describe the gas 
kinematics close to the (proto)star. 



4.2.2. 6.7 GHz Methanol Maser 



Figure [3] (middle panel) shows the total-power spectrum of the 6.7 GHz methanol 
masers observed on 2008 March 18 toward G28.87+0.07 using the Effelsberg 100-m 
telescope. The 6.7 GHz methanol maser spectrum, with two narrow features at Vlsr of 



abo ut 97 and 105 
that Caswell et al. 



t m s \ resembles that of the water masers (Figure [31 upper panel). Note 
(1995|) detected only the emission at Vlsr of 105 km s -1 with the Parkes 
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64-m telescope. 

We imaged the range of LSR velocities from 95 to 107 km s _1 with the EVN and 
detected two distinct methanol maser features, both with Vlsr close to 105 km s -1 . Feature 
#1 is detected at the first, third, and fourth epoch, whereas feature #2 is detected at 
the third and fourth epoch only. Due to adverse atmospherical conditions and telescope 
failures, no 6.7 GHz maser signal was detected at the second epoch. Individual feature 
properties are presented in Table HI Positional offsets are relative to feature #1, and the 
absolute positions of feature #1 at the three epochs of detection are given in Table [2j The 
two methanol maser features lie at the western border of the distribution of the red-shifted 
water masers, ~0'/2 to the southwest of the VLA 3.6 cm continuum peak (see Figure H]). 
Inspecting our images, we could not find any reliable emission at the velocity of ~97 km 
s _1 , in correspondence of the weakest feature detected in the total-power spectrum. The 
mean rms noise on the channel maps was 20 mJy beam -1 , and we could have missed that 
emission if falling below the 5a detection threshold. Because of the narrow line width 
(<1 km s -1 ) and the offset in Vlsr from the systemic velocity (101.0 km s _1 ) of the HMC, 
we tend to exclude that the weak methanol feature at Vlsr ~97 km s -1 is thermal emission 
from the HMC resolved out by the EVN observations. 

After correcting for the apparent motion, the absolute proper motion for feature #1 is: 
V a = —27.8 ± 1.0, Vg = —45.7 ± 4.6 km s -1 . As discussed in Sect. I4.2.1[ we have used this 
vector to correct the absolute proper motions of the water masers. The relative motion of 
the methanol maser feature #1 with respect to the "center of motion" of the water masers 
is found to be: v a = 1.7 ± 6.7, vg = —15.6 ± 8.0 km s _1 . This is consistent with the results 
from other VLBI observations which show that 6.7 GHz masers move significantly slower 
than water masers, and justifies employing the methanol maser absolute motion to correct 
the water maser absolute motions. 
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4-2.3. 1. 665 GHz OH Maser 



Forster fc Caswelll (119991 ) observed the HMSFR G28.87+0.07 with the VLA A-B 
hybrid configuration (HPBW of 5'.'5 x 1'.'2) in the left circular polarization (LCP) band 
and detected eight 1.665 GHz OH maser spots. Their brightest spot had a flux density of 
9.85 Jy at the V LSR of 102.8 km s -1 . 

Figure [3] (lower panel) shows the T sys -weighted mean of the VLBA antenna total-power 
spectra of the 1.665 GHz OH maser emission in the LCP and RCP bands. Three strong, 
narrow features are visible in the LCP band over the LSR velocity range from 102 to 108 km 
s -1 . Two weak features are detected in the RCP band, one of them a t an LSR velocity 
of ~1 08 km s -1 . These findings are similar to the results obtained by 



Caswell & Haynes 



( 119831 ) with the Parkes 64-m telescope. 



We imaged the range of LSR velocities from 102 to 110 km s -1 , and identified 3 distinct 
OH maser features in the LCP band. In the RCP band, by imaging maser data after direct 
phase-referencing to the calibrator J1834— 0301, no signal was detected above a 5<r level 
of 0.1 Jy beam -1 . Individual feature properties are presented in Table El The positional 
offsets are relative to feature #1, whose absolute position is give n in Table El The derived 



Forster k Caswelll (11999h . 



absolute position is offset by ~0'/9 from the reference position of 
Maser intensities range from 0.7 to 1.9 Jy beam -1 . The LSR velocities vary from 102.8 km 
s -1 for the most blue-shifted feature (#1), to 105.5 km s -1 for the most red-shifted feature 
(#3). The FWHM line width of individual maser features ranges from 0.4 to 1.4 km s -1 . 
Positions and LSR velocities of hydroxyl maser features are plotted in Figure HI The three 
detected features lie ~ 0'/6 to the southeast of the VLA 3.6 cm continuum source and are 
distributed along a line extended ~ 0'/3. 
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4.3. 24.5 /im Infrared Emission 

Figure [7] shows an image of the 24.5 /im continuum emission obtained with 
Subaru/COMICS at an angular resolution of ft' 75. One can see that the structure of the 
emitting region is quite complex. Beside two components associated with the radio sources 
"HMC" and "A" with a total flux of ~80 Jy, another fainter source (~3 Jy) is visible to the 
south, offset by ~4" from "A" . 

The emission from the northern region is clearly resolved and appears to be roughly 
associated with the two radio objects, although one cannot rule out the existence of 
additional, fainter unresolved sources contributing to the extended halo around the two 
main sources. In particular, the 24.5 /im emission from "HMC" seems to split into two 
sub- components in the E-W direction, with comparable flux densities: ~32 Jy for the 
eastern and ~43 Jy for the western component. These could trace two distinct objects or 
diffuse emission surrounding a single source at the origin of the "HMC" radio continuum. A 
third, more intriguing possibility is that one is looking at an "hour-glass" shaped nebulosity 
centered on the embedded source at the origin of the free-free emission. 

In the latter scenario, the source lies at the center of a flattened core lying approximately 
along the line of sight. This dusty core absorbs most of the IR photons emitted at 24 /im 
towards the observer, while the IR radiation can escape along the minor axis of the core, 
lying close to the plane of the sky. We have estimated an optical depth of ~ 4 for the HM C 



at 24 /im, assuming the H2 colum n density of 4.2 x 10 23 cm 2 obtained by 



( 120081 ) and the opacity adopted by 



Whitney et al. 



Furuya et al 



( I2003a| ) - see their Figure 1. Such a value 



could be sufficiently large to explain the lack of emission towards the embedded star and 
determine the "hour-glass" shape of the IR emission. Interestingly, the E-W orientation of 
this putative bipolar nebula is also roughly consistent with the direction of the H2O maser 
jet in Figure EJ suggesting that one could be observing a bipolar flow piercing through a 
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dense core, and revealed on different scales in different tracers. However, such a scenario, 
albeit appealing, is not consistent with the inclination of the maser jet with respect to the 
plane of the sky. In fact, the blue-shifted emission lies to the E (see Figure E]) and this 
implies that the eastern lobe should be pointing towards the observer and thus be brighter 
than the western lobe in the IR. Looking at Figure one can see that this is not the case, 
being the emission slightly stronger to the west (see above). 

In conclusion, we consider both the "hour-glass" nebula and the double source 
hypotheses too speculative on the basis of the observational evidence collected so far, and 
in the following we make the conservative assumption that the IR emission is associated 
with "HMC", with a flux of -75 Jy. 

The southern source "A" has a 24.5 fim flux density of —5 Jy. Also in this case, the 
association between the IR emission and the radio source is questionable, because the two 
are offset each other by ~1". However, this is the astrometrical error of the IR image and 
we thus believe that such an offset is not sufficient evidence against a common origin for 
the IR and radio emission toward "A". Therefore, we conclude that both emissions can 
originate from the same object. 



4.4. Spectral Energy Distribution 



Figure M shows the SED of G28.87+0.07 from 3.6 /im to 1.1 mm. Beside our COMICS 



data, we have used data from the f ollowing surveys: Ga 



Survey Extraordinaire 



MSX flEgan et al. 



(GLIMPSE; 



2003J) . IRAS flNeuge 



Benjamin et a 



Plane Survey (Hi-GAL; iMolinari et al. 



of the Galaxy (ATLASGAL; 



bauer 



et al. 



2010a 



Schuller et al. 



actic Legacy 



2003), MIPSGAL 



nfrared Mid-Plane 



Carey et al. 



2009|), 



19841 ). the Herschel infrared Galactic 



bj), the APEX Telescope Large Area Survey 



20091 ). and the Bolocam Galactic Plane Survey 
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(BGPS; iDrosback et al.ll2008l ). For the Herschel/Hi-GAL data, the most recent release 
(obtained in July 2011 with HIPE version 7.0.0) w as used, resulting from data reduction 



201 ll ). where the most relevant image artifacts 



with the RomaGal software ( ITraficante et al 
have been removed. 

Apart from the IRAS and MSX flux densities, which have been taken from the 
corresponding Point Source Catalogues (sources IRAS 18411-0338 and G028.8621+00.0657), 
all the other values have been calculated by integrating the emission inside the dashed 
contour in Figure [9j The latter was chosen in such a way to include the whole emitting 
region at all wavelengths. In this way we overcome the problem of different angular 
resolutions at different wavelengths. In the calculation of the flux the background emission 
was subtracted assuming that the background flux inside the contour is in all points equal 
to the typical flux density measured along the contour itself. The flux densities used for the 
SED are given in Table [6j 

Figure [9] illustrates the pc-scale structure of the G28.87+0.07 star forming region, 
imaged with Herschel from 70 to 500 /im and with APEX at 870 /im, and demonstrates 
that from the mid-IR to the sub-mm the continuum emission is dominated by a compact 
source surrounded by a weaker halo. 



We hav e fitted the SED with th e radiative transfer model developed by 



(120070 and 



Whitney et al. 



Robitaille et al. 



( l2003al Jbl) (using the SED fitting tool available on 



http://caravan.astro.wisc.edu/protostars/), which assumes a ZAMS star with a cir- 
cumstellar disk, embedded in an infalling envelope, and allows derivation of a number of 
physical parameters. While we are aware that this model cannot adequately reproduce a 
complex situation such as that of the multiple MYSOs in G28.87+0.07, we believe that 
the fit is sufficiently good to guarantee good estimates of integrated quantities such as the 
luminosity and the mass of the envelope. The best-fit values are, respectively, 2.3 xlO 5 L Q , 
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and 3.0 xlO 3 M & . Note that in the fitting procedure, the IRAS fluxes were considered upper 
limits due to the limited angular resolution, whereas the MIPSGAL 24 fim flux is in fact a 
lower limit because the corresponding image is saturated and only the non saturated pixels 
have been considered in the calculation. Also, the distance was kept fixed to the nominal 
value of 7.4 kpc. The best-fit value for the visual foreground extinction is 46 mag. 



5. DISCUSSION 

5.1. Nature of the Radio Continuum 

In this section, we focus on the cm-continuum source labeled "HMC" in Table HJ found 
in good correspondence with the positions of the 22 GHz H2O, 6.7 GHz CH3OH, and 
1.665 GHz OH masers in G28.87+0.07 (see Figure Hj). The continuum spectrum (shown in 
Figure [2]) is consistent with free-free emission from an ultracompact (UC) H II region, with 
the turnover frequency in the range 10-20 GHz, and ionized by a Lyman continuu m rate 



of 7.9 xlO 45 s l . This would correspond to a ZAMS star of spectral type B1-B0.5 (jPanagia 



19731 ). with a luminosity of L sta r ~ 8 x 10 3 L & , and a mass of M star ~ 10 M & . It is worth 



mentioning that these values would b e different if mo re recent studies than the seminal 



article by Panagia were adopted (e.g., 



Crowther 



20051 ). However, the basic result that one 
is dealing with an early B-type star with a luminosity of ~ 10 4 L & would remain valid. 

We note that the fit to the continuum spectrum requires a very small size of the 
UC H II region, of ~0'/04 or ~300 AU, which appears in contradiction with the fact that 
our A-array 1.3 cm observations partly resolve the continuum emission (see Tabled]). 

Since we considered a homogeneous distribution of ionized gas in our fit, one possibility 
is that the H II region is not homogeneous. The southern spur identified in the A-array 
3.6 cm map of the "HMC" source (see Figure H]), could indeed suggest the existence of a 
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density gradient toward the south. One may speculate that being source "A" more extended 
and perhaps more evolved, the gas in its surroundings is less dense. In this scenario, the 
north-south streamlines of water masers on either sides of the southern spur, might trace 
shocks at the interface between the ionized gas of "HMC" (expanding to the south) and the 
surrounding molecular gas. 

We also note that the continuum spectrum of "HMC" between 1.3 and 6 cm can be 
also fitted with a power law S v oc v a , w ith a = 0.6 (see Figure [2]). This spectral index is 



typical of thermal jets ([Reynolds 



19861 ) and provides us with an alternative explanation 
for the origin of the free-free emission. However, a caveat is in order. As discussed in 
Section H~Tl the slope of the spectrum between 3.6 cm and 1.3 cm could be significantly 
affected by the inteferometer resolving out part of the emission at the shorter wavelength. 
If the spectrum is steeper than estimated from our data, the slope could be still consistent 
with a thermal jet , but also with that of an H II region with a density gradient (see e.g., 



Franco et al 



2000 



We conclude that it is impossible to decide the nature of the free-free emission on 
the basis of the information collected so far. Notwithstanding this caveat, as a working 
hypothesis we prefer the jet interpretation, because this permits to interpret the continuum 
and maser emission in the same context, as discussed in the following section. 



5.2. The Water Maser Flow 

In this section we propose that both the continuum and the water maser emission are 
associated with a thermal jet impinging against the surrounding quiescent molecular gas. 
This idea is based on two observational findings. First of all, Figure |5] shows that, at larger 
distances from the center of the water maser distribution, maser (relative) velocities appear 
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to be larger and better collimated along the northeast-southwest direction. Secondly, the 
axis of the powerful, molecular outflo w mapped at arcsec ond scales toward G28.87+0.07 



has a northeast-southwest orientation (IFuruya et al. 



20081 . Figure 5). This it is plausible to 



assume that the two, fast moving, water maser clusters "Jb" and "J r ", located, respectively, 
at the northeast and southwest edge of the maser distribution, are excited by the same jet 
responsible also for driving the large-scale molecular outflow. One could object against the 
"jet" interpretation, that the whole pattern of water maser velocities cannot be explained 
only in terms of a collimated outflow. In fact, at smaller distances from the center of the 
maser distribution, water maser velocities point to the north or south directions, and form 
large angles with the direction of the putative jet. We propose that the scattered pattern of 
proper motions presented by the red-shifted water masers can result from the interaction of 
the jet with dense circumstellar material, which deflects and increase the opening angle of 
the flow velocities. 

If the free-free continuum source "HMC" is a jet, the position of the (proto)star 
powering it may be offset from the free-free continuum peak. If a jet powered by the 
(proto)star drives the water maser motions, and, in particular, the expansion to the south 
and west of the narrow strip of red-shifted water masers, the (proto)star could be located 
at an intermediate position between "HMC" and the red-shifted maser strip. This point is 
not far (within tens of mas) from the position of the center of motion of the water masers 
(cross in Figure |5]). Considering the large uncertainties in deriving the center of water 
maser expansion, in the following discussion we can safely take the center of motion as a 
reasonable estimate of the (proto)star location. 

In the jet scenario, one can estimate the momentum rate of the jet from the 
average distance of the H2O masers from the MYSO and the average maser velocity. 



Assuming a H 2 pre-shock density of 10 8 cm 3 , as predicted by models (jElitzur et al 
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19891 ). the momentum rate in the water maser jet is given by the expression P = 
1.5 x 10~ 3 V 2 Rl 00 (Q / An) M Q yr -1 km s _1 . Here Vio is the average maser velocity in units 
of 10 km s _1 , i?ioo the average distance of water masers from the YSO in units of 100 AU, 
and Q the solid angle of the jet. This expression has been obtained by multiplying the 
momentum rate per unit surface transferred to the ambient gas (riH 2 m H 2 V 2 ), by QR 2 , 
under the assumption that the jet is emitted from a source at a distance R from the masers 
within a beaming angle Q. To estimate the average distance and velocity of the water 
masers tracing the jet in G28.87+0.07, we use all the maser features belonging to the two 
fast moving clusters "J&" and "J r " (see Figure E]). The resulting average distance is 122 mas 
(or 902 AU) and the average speed is 43 km s _1 . Using these values, the jet momentum 
rate is P = 2.2 (to /Air) M Q yr" 1 km s _1 . 

For shock induced ionization in a thermal jet, the momentum rate of the jet can 



be estimated 



Anglada et al. 



rom t 



re continuum flux assuming the emission is optically thin (see 
(|1996[ )). One has: F v d 2 = 10 3 ' 5 (f2/47r)P, where F v is the measured 
continuum flux density in mJy, P the jet momentum rate in M & yr~ x km s -1 , Q the jet 
solid angle in sterad, and d the source distance in kpc. Using the flux of 1.2 mJy measured 
at 1.3 cm with the VLA C-array for source "HMC", and the distance of 7.4 kpc, one derives 
P = 2 x 10~ 2 {VL/Ati)- 1 M q jt- 1 km s" 1 . 

Note that the momentum rates estimated from the maser motion and the continuum 
flux have a different dependency on the jet solid angle, which can thus be estimated by 
equating the two momentum rate estimates. 

We obtain an expression for the jet solid angle: 2.2 (Q/Att) = 2 x 10~ 2 (Q/Att)^ 1 . 
We find Q = 1.2 sterad, i.e. a jet semi-opening angle 6j = 35°. Correspondingly, the jet 
momentum rate is P = 0.2 M & yr _1 km s _1 . This value is lar ge enough to j ustify the 
momentum rate of several 10~ 2 M & yr _1 km s _1 , estimated by 



Furuya et al. 



O2008|) for 
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the large-scale molecular outflow. This result supports our hypothesis that the small-scale 
thermal and H.?0 maser je t is su fficiently powerful to feed the large-scale molecular outflow 



observed by 



Furuya et al. 



(120081 ). We note that the derived parameters for the jet in 



G28.87+0.07 are comparable to those of the few collimated outflows o bserved up to now 



in massive star-formin g regions. For compar ison, the jets in HH 80-81 (M arti et al 



in IRAS 16547-4247 (IRodrfguez et al. 



1998) 



20081 ). and in IRAS 20126+4104 (Moscadelli et al. 



201 ll ) have momentum rates of the order of 10 1 M Q yr 1 km s 1 and opening angles of 



20-30°. 



5.3. Methanol Maser Environment 



The two detected 6.7 GHz methanol maser features are approximately aligned along 
the same northeast-southwest line defined by the two maser jet clusters "J&" and "J r " , the 
peak of the VLA A-array 1.3 cm, and the center of motion (see Figure H]). Figure E] shows 
that the total-power spectra of the 6.7 GHz methanol and 22 GHz water masers are similar, 
both presenting a strong red-shifted and a much weaker blue-shifted component. The 
positions of the methanol maser features roughly along the jet axis, and the similarities in 
the spectral shape of the water and methanol maser emissions, make us speculate that also 
the 6.7 GHz methanol masers could originate in the (proto) stellar jet, albeit associated with 
more quiescent gas than that traced by the wate r masers. Maser VLB! observations toward 



2011 



Sanna et al. 



2010bh 



the sources IRAS 20126+4104 and G23.01-0.41 (iMoscadelli et al. 
have also indicated that the 6.7 GHz masers might trace entrained gas at the interface with 
a jet. 

Figure H] clearly shows that the blue-shifted lobe of the water maser jet is sampled 
much more sparsely, as it presents much less features than the red- shifted lobe. This is 
consistent with the blue-shifted line being significantly weaker than the red-shifted one 



-28 - 



in the total power spectrum (see Figure |3l upper panel). This fact suggests that the 
circumstellar material is not isotropically distributed around the (proto)star, but it is 
mostly concentrated on the side of the red-shifted lobe. That would also account for the 
stronger emission line in the total-power spectrum of the CH3OH masers at the red-shifted 
velocities. 



5.4. OH Maser Environment 



The spectrum of the 1.665 GHz maser is dominated by three narrow LCP features (see 
Figure |3l lo wer panel) , and is very sim ilar to the spectrum of the 1.667 GHz maser transition 



observed by 



Caswell fc Hayned (119831 ). These authors find that most OH interstellar masers 



are highly polarized, either RCP or LCP, with the same predominant polarization for both 
the 1.665 GHz and 1.667 GHz transition. They propose an explanation based on large-scale 
velocity and/or magnetic field gradients across the masing region. 

Figure [1] shows that the 1.665 GHz OH maser features are found significantly offset 
(« &!6) from the compact VLA source "HMC" and the center of water maser activity, and 
are closer to the more extended source "A". A possible interpretation is that OH masers 
trace more external, less dense parts of the molecular core surrounding the (proto)star 
responsible for the excitation of the water and methanol masers. Observations indicate 
that the density of massive proto-stellar c ores decreases from t he center outward, with a 
power-law index of about 1.6 (n oc r~ 16 ) (iBeuther et al.ll2002ai ). In G28.87+0.07, the ratio 
of the average distance (from the "center of motion" of water masers) of the water (« O'.'l) 
and OH (« 0'.'6) masers, is about 6, implying a decrease in gas density by a factor of 
about 20. That agrees with maser ex citation models, which predict a H? pre-shock density 



id 



of 10 6 -10 7 cm" 3 for OH 1.665 GH 
10 8 cm~ 3 for H 2 22 GHz masers (jElitzur et al 



Cesaroni fc Walmslev 



1991 



Cragg et al. 



2002), and 



19891 ) . respectively. 
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5.5. The Nature of Source "HMC" 

So far, we have explored the properties of source "HMC" making use of the maser, 
radio, and IR continuum emission. Now, we want to use our findings to shed light on the 
nature of the MYSO(s) associated with this source. 

We have seen that the radio continuum emission and the water maser features suggest 
the presence of a jet with a large momentum rate, on the order of ~ 10 _1 M & km s _1 yr _1 . 
Values that high are consistent with the jet being powered by a high-mass (proto)stars, as 
indicated, e. g., by the relat i onship between outflow momentum rate and MYSO luminosity 



obtained by iBeuther et al 



(l2002bf ). In particular, from their Figure 4b, one can conclude 
that the MYSO must have a luminosity of at least several 10 3 L & . Alternatively, if the 
radio emission were originating from an H II region, we have seen (Section I5.ip that the 
ionizing star should have a luminosity < 10 4 L & . Clearly, both the thermal jet and the H II 
region hypotheses imply similar luminosity for the associated MYSO. The question we want 
to investigate here is whether such a luminosity is consistent with the results obtained from 
the IR data. 

In Section 14.41 we have found that the total luminosity of the region is on the order of 
2 x 10 5 L & . However, this includes all MYSOs in a radius >10" around "HMC", whereas 
we need to estimate the contribution of the sole "HMC" , to achieve a consistent comparison 
with the luminosity estimate obtained from the water maser and/or radio continuum. Our 
24.5 fim continuum image can be used for this purpose. As shown in Figure [7], the "HMC" 
source is dominating the emission at this wavelength, contributing to ~90% of the total 
flux density, and a similar conclusion holds also for the GLIMPSE images. Due to the lack 
of sub-arcsecond imaging at longer wavelengths, it is difficult to establish whether also 
the flux density measured close to the peak of the SED is mostly due to source "HMC". 
Nonetheless, the Hi-GAL/Herschel data suggest that the 70 /im emission peaks closer to 
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source "HMC" than source "A" . A word of caution is in order here. The angular resolution 
of the Hi-GAL data at 70 /iin is estimated to be about 10" x 9" (after reduction with the 
Roma-Gal software), while the astrometrical error is ~3'.'5 (the latter results from aligning 
the PACS 70 /zm images to the MIPSGAL 24 /im images by matching suitable point sources 
detected at both wavelengths). Although large, these values still permit to associate the 
70 fim. emission with source "HMC". In fact, the angular separation between the 3.6 cm 
continuum source "HMC" and the 70 /j,m peak is 3'.'7, consistent with the positional error, 
whereas source "A" lies 6'/5 from the 70 /im peak. This can be appreciated in Figure [TOj 
where one sees that the 70 /im emission peaks to the NW of "HMC". Therefore, in the 
following we assume that the bolometric luminosity estimated from the SED is likely 
dominated by the "HMC" source and conclude that the MYSO(s) embedded in the HMC 
may be as luminous as ~ 2 x 10 5 L & . 

Notwithstanding the large uncertainties of all our luminosity estimates, the value 
obtained from the IR images exceeds by an order of magnitude that derived from the maser 
and radio continuum data. This inconsistency can be explained if "HMC" is hosting a 
multiple stellar system rather than a single star. One can get an idea of how this may 
affect our estimate of the stellar properties (ma ss, etc.) by assuming that the luminosity of 



2 x 10 5 L & is emitted by a stellar cluster with a 



Miller fc Scald ( 119791 ) initial mass function. 



One finds that the luminosity of the most massive star may be as small as ~ 3 x 10 4 L , 
depending on the statistical distribution of the stellar masses in the cluster. Given the 
uncertainties, this value is only marginally greater than the luminosity estimate from the 
radio continuum and maser data. 

The same method provides us with the total mass of the stellar cluster, which ranges 
from 5 x 10 2 M & to 2 x 10 3 M , depending on the statistical distribution of the stars in 
the cluster. This stellar mass can be divided by the total mass of the associated molecular 
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clump, 3 x 10 3 M & (see Section E3]) , to calculate the star formation efficiency, 17%-67%. 



Albeit large, values like these are not implausible, as we are considering the part o: 



molecular cloud where star formation concentrates (ILada fc Lada 



2003 



Bonnell et al. 



the 



201lh . 



In conclusion, we believe that the observational evidence collected so far for source 
"HMC" suggests that one is observing a young star of 20-30 M , still in the main accretion 
phase and surrounded by a rich stellar cluster. Such a conclusion holds independently of 
the scenario adopted to explain the free- free emission (thermal jet or weak H II region) 
and raises a question: why only faint free-free emission is detected from this object? A 
~ 3 x 10 4 L & star is expected to ionize the surrounding gas emitting a flux density of ~2 mJy 
(in the optically thin limit), whereas towards source "HMC" less than 1 mJy is measured. 



A possible explanation is that the H II region is at least 



part 



Ketol (12002 . 



quenched by accretion 



2003 !) . Although n o 



2008|), 



from the parental core, as depicted, e.g., in the model by 

evidence of accretion is seen in a typical infall tracer such as HCO + (iFuruya et al. 
this fact is inconclusive because the majority of the detected HCO + (1-0) might be part of 
the outflowing gas, thus hindering the detection of infall. 

Finally, it is also worth mentioning that according to recent theoretical models 



(IHosokawa et al. 



20101 ). large accretion rates are bound to bloat the protostellar radius and 
decrease the corresponding temperature and Lyman continuum flux. This could be an 
alternative explanation of the lack of a (strong) H II region. 



6. SUMMARY 



We observed the HMSFR G28.87+0.07 with the VLBI phase-referencing technique 
in three powerful maser transitions: 22.2 GHz H2O, 6.7 GHz CH3OH, and 1.665 GHz 
OH. In addition, we also performed VLA observations of the radio continuum emission at 
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1.3 and 3.6 cm with both the A- and C-array, and Subaru/COMICS observations of the 
mid- infrared continuum emission. We also made use of data from the Hi-GAL/Herschel, 
ATLASGAL/APEX, as well as other surveys of the Galactic plane. From all these 
observations, we draw the following conclusions: 

1. The bipolar distribution of l.o.s. velocities and the general pattern of observed relative 
proper motions indicate that the water masers are tracing expansion with an average 
velocity of ~24 km s _1 from a point on the sky close to the maser "center of motion" . 
While at larger distance from the "center of motion" , maser velocities are higher and 
better collimated along the northeast-southwest direction, at closer distances water 
masers show smaller speeds and a larger scatter in the velocity orientations. We 
interpret these facts as the result of the interaction of a (proto)stellar jet with dense 
circumstellar material. 

2. Two continuum sources (named "HMC" and "A"), separated by ~3'.'l along the 
north-south direction, are observed with the VLA at 3.6 and 1.3 cm. The "HMC" 
source is compact and spatially associated with the observed H 2 0, CH 3 OH and OH 
masers. Whether the free-free emission arising from this source is due to an H II 
region or a thermal ionized jet cannot be unambiguously established with the present 
data. However, we believe that the jet hypothesis is more consistent with the direction 
of the water maser motions, which outline expansion in a bipolar flow. Moreover, the 
orientation and momentum rate of the putative jet are compatible with those of the 
molecular outflow detected at arcsecond scales. We hence argue that the jet traced by 
the water masers is also powering the large-scale outflow. 

3. The combined information obtained from sub-arcsecond imaging with Sub- 
aru/COMICS at 24.5 /irn and lower angular resolution images from the Hi- 
GAL/Herschel and ATLASGAL/APEX surveys, permit to estimate with good 
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accuracy the luminosity (2 x 10 5 L Q ) and gas mass (3 x 10 3 M ) of the star forming 
region of G28.87+0.07, and establish that ~90% of the luminosity is coming from the 
radio source "HMC". We conclude that this source must contain multiple stars, with 
the most massive being at least as luminous as 3 x 10 4 L Q . The lack of an associated 
Hn region indicates that the 20-30 M star must be still undergoing heavy accretion 
from the surrounding envelope. 

We thank the anonymous referee for careful reading of the text and the constructive 
criticisms that improved the presentation and analysis of our results. This work was 
supported by the Chinese NSF through grants NSF 11133008, NSF 11073054, NSF 
10733030, and NSF 10703010. 
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Fig. 1. — Maps (gray scale) of the VLA continuum emission toward G28.87+0.07 at 1.3 and 
3.6 cm. (a) 3.6 cm emission with the VLA C-array. (b) 1.3 cm emission with the VLA 
C-array (c) 3.6 cm emission with the VLA A-array (d) 1.3 cm emission with the VLA A- 
array. Squares and circles report the positions of the 1.665 GHz OH and the 6.7 GHz CH 3 OH 
masers, respectively, detected in this paper (see Section. 14.2.21 and I4.2.3H . In each panel, the 
contour levels range from 30% to 90% of the peak emission at multiples of 10%. The restoring 
beam is shown in the lower, left or right corner of each panel. A rectangle overlaid on the 
image indicates the field of view expanded in the adjacent panel. The associated labels 
indicate the names of the continuum sources identified by us (see Section fl~Tl and Tabled]). 
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Fig. 2. — Spectral energy distribution of the VLA component "HMC" toward the HMSFR 
G28.87+0.07. Dots and error bars report the values and the associated errors (20%) of 
the fluxes at 1.3 and 3.6 cm measured with the VLA C-array (see Tabled]). The triangle 



indic ates the upper limit of the flux at 6 cm obtained with the VLA B-array (IPurcell et al. 



20081 ). The linear fit to these measurements, indicated by the dotted line, produces a spectral 
index a = 0.6. The solid line shows the best fit of the measured fluxes with a model of a 
homogeneous H II region, yielding a Lyman continuum rate of 7.9 x 10 45 s _1 and a H II region 



radius of 0"04. 
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Fie. O. Total-power spectra of the H2O, CH3OH, and OH masers toward G28.87+0.07. Upper panel: system-temperature 
(T sys ) weighted average of the 22.2 GHz total-power spectra of the 9 VLBA antennas observing on 2006 April 23. Middle 
panel: Effelsberg total-power spectrum of the 6.7 GHz methanol maser emission on 2008 March 18. Lower panel: T S y S -weighted 
average of the 1.665 GHz total-power spectra of the 9 VLBA antennas observing on 2007 April 28. Continuous (black) and 
dashed (red) lines are used to distinguish between right (RCP) and left (LCP) circular polarizations, respectively. The dotted 
line crossing the spectra denotes the systemic velocity (V S ys) inferred from NH3 measurements. 
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Fig. 4. Spatial distribution of 22 GHz H 2 (triangles), 6.7 GHz CH 3 OH (dots), and 1.665 GHz OH (squares) masers in 

G28.87+0.07. Different colors are used to indicate the maser LSR velocities, according to the color scale on the righthand side 
of the plot. Upper Panel: global view of the region where the three maser species distribute. The solid contours give the 3.6 cm 
continuum emission observed with the VLA A-array. Contour levels range from 20% to 90% of the peak emission (0.35 mjy 
beam -1 ) at multiples of 10%. Dashed rectangles enclose the two north-south streamlines of water masers observed on either 
side of the VLA 3.6 cm continuum. A dotted rectangle overlaid on the image indicates the field of view expanded in the right 
panel. Lower Panel: zoom on the area where most of water masers concentrate. The solid contours give the 1.3 cm continuum 
emission observed with the VLA A-array. Contour levels range from 40% to 90% of the peak emission (0.62 mjy beam" 1 ) at 
multiples of 10%. 



-42 - 




46*220 



RIGHT ASCENSION (J2000) 




t=3 




46"226 



46fe24 



km s 
93 
94 
95.1 
96.1 
97.1 
98.2 
99.2 
100.2 
101.2 
102.3 
104.3 
105.3 
106.4 
107.4 
10B.4 
109.4 
110.5 
111.5 
112.6 



RIGHT ASCENSION (J2000) 

Fie. 5. Proper motions of 22 GHz H2O masers (triangles) relative to their "center of motion", and spatial distribution of 

6.7 GHz CH3OH masers (dots) in G28.87+0.07. Different colors are used to indicate the maser LSR velocities, according to the 
color scale on the righthand side of the plot. The cones indicate the 3-D velocities of the water maser features relative to the 
"center of motion" (marked with a cross). Points without an associated cone have been detected only over one or two epochs and 
the associated proper motion cannot be computed or is considered unreliable. The cone opening angle gives the la uncertainty 
on the proper motion direction. The length of the cone is proportional to the velocity, with the amplitude scale indicated in the 
lower left corner of the figure. Open dotted rectangles mark two clusters of water masers located at the northeast and southwest 
edge of the distribution (labeled "Jj," and "J r ", respectively), moving fast and close to the northeast-southwest direction. The 
VLA 1.3 cm continuum emission is plotted with solid contours. Contour levels range from 40% to 90% of the peak emission 
(0.62 mjy beam -1 ) at multiples of 10%. Upper Panel: view of the whole distribution of water maser proper motions. Lower 
Panel: zoom on the most crowded area of water maser concentration, to the west of the VLA 1.3 cm continuum. 
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Fig. 6 . Absolute proper motions of water masers in G28.87+0.07, derived by applying different corrections. Upper 

Panel: Correcting for the apparent motion due to the combination of the parallax, the solar motion with respect to the LSR, 
and the differential Galactic rotation (adopting the flat rotation model indicated in Section |3.2| |. Lower Panel: Correcting for 
the absolute proper motion of feature #1 of the CH3OH masers. Symbols, colors, and contour levels have the same meaning 
as in Figure [5] The cones indicate the 3-D velocities of the water maser features. Points without an associated cone have been 
detected only over one or two epochs and the associated proper motion cannot be computed or is considered unreliable. The 
cone opening angle gives the la uncertainty on the proper motion direction. The length of the cone is proportional to the 
velocity, with the amplitude scale indicated in the lower left corner of the figure. 
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Fig. 7. — Continuum emission maps toward G28.87+0.07 at centimeter and infrared wave- 
lengths. The grey image in the background is the 24.5 /im emission from Subaru/COMICS. 
The white and black contours indicate the centimeter continuum emission at 3.6 and 1.3 cm, 
respectively. The values of the contours are the same as in Figure [TJ 
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Fig. 8. — Spectral energy distribution of the VLA component "HMC" toward G28.87+0.07 
from 3.6 /im to 1.1 mm. Dots indicate real values, whereas triangles with the vertex up 
and down denote lower and upper limits, respectively. Section 14.41 lists the infrared tele- 
scopes used to reconstr uct the SEP. W e have fitted the SED with the radiative trans- 
fer model developed by 



Robitaille et al. 



(120071 ) (using the SED fitting tool available on 



http:/ /caravan. astro. wisc.edu/protostars/ ) 
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Fig. 9. — Each panel shows the pc-scale structure of the G28.87+0.07 star-forming region 
at a different infrared wavelength (indicated in the upper left corner of the panel). Images 
from 70 to 500 fim are observed with Herschel, and the map at 870 /zm is taken with APEX. 
The color scale corresponds to increasing intensity from dark purple to white, as indicated 
in the wedge to the right. The observing beam is reported in the lower right corner of each 
panel. The dashed white circle denotes the area of integration to calculate the G28.87+0.07 
flux at different infrared wavelengths. 
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Fig. 10. — Overlay of the 3.6 cm continuum map (contours; same as in Figure Cj) on the 
Hi-GAL/Herschel image at 70 fim. Note how the far-IR emission peaks close to the free-free 
source "HMC". 
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Table 1. Radio continuum emission in G28.87+0.07. 



Peak position 



Label 


Array 


A 


HPBW 


PA 


Image rms 


R.A.(J2000) 


Dec.(J2000) 


Fpcak 


F int 






(cm) 


(" x ") 


(°) 


(mjy beam -1 ) 


m s \ 


(° ' ") 


(mjy beam -1 ) 


(mjy) 


HMC 


VLA-A 


1.3 


0.122 x 0.093 


-17 


0.06 


18 34 46.23 


-03 35 29.8 


0.62 


0.64 




VLA-C 


1.3 


1.053 x 0.755 


-9 


0.05 


18 43 46.23 


-03 35 29.8 


1.03 


1.20 




VLA-A 


3.6 


0.33 x 0.24 


-1 


0.03 


18 43 46.23 


-03 35 29.9 


0.35 


0.42 




VLA-C 


3.6 


2.62 x 2.11 


5 


0.02 


18 43 46.25 


-03 35 29.8 


0.53 


0.55 


A 


VLA-A 


1.3 


0.122 x 0.093 


-17 


0.06 






<0.18(3ct) 






VLA-C 


1.3 


1.053 x 0.755 


-9 


0.05 


18 43 46.30 


-03 35 32.8 


0.39 


0.60 




VLA-A 


3.6 


0.33 x 0.24 


-1 


0.03 






<0.09(3ct) 






VLA-C 


3.6 


2.62 x 2.11 


5 


0.02 


18 43 46.29 


-03 35 32.7 


0.81 


0.98 



Note. — Column 1 gives the source label; columns 2 and 3 report the VLA array configuration and the observing band, respectively; 
columns 4 and 5 the (major and minor) HPBW size and the PA of the observing beam, respectively; column 6 the image rms; columns 7, 
8, 9, and 10 the absolute position, the peak intensity, and the integral intensity of the emission peak, respectively. The PA of the 
beam is defined East of North. 
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Table 2. Positions and Brightness. 



Maser 


Feature 


epoch 


R.A. (J2000) 


Dec. (J2000) 


Fpcak 


NW beam 








m s) 


(° ' ") 


(Jy beam ) 


(mas X mas 


@ deg) 


H 2 


1 


2006 Apr 23 


18 43 46.22467 


-03 35 29.8161 


57.4 


1.2x0.4 @ 


-17 






2006 Jun 30 


18 43 46.22465 


-03 35 29.8173 


50.3 


1.3x0.4 @ 


-16 






2006 Sep 28 


18 43 46.22455 


-03 35 29.8190 


32.0 


1.4x0.5 ® 


-3 






2007 Jan 18 


18 43 46.22448 


-03 35 29.8212 


20.7 


1.5x1.0 @ 


20 


CH3OH 


1 


2006 Feb 28 

2007 Mar 18 


18 43 46.22079 


-03 35 29.8379 


0.044 


9.9x5.6 @ 


5 






2008 Mar 18 


18 43 46.22031 


-03 35 29.8521 


0.246 


12.7x5.2 € 


j 45 






2009 Mar 17 


18 43 46.22011 


-03 35 29.8595 


0.175 


12.6x5.9 t 


5 13 


OH 


1 


2007 Apr 28 


18 43 46.26771 


-03 35 30.4697 


1.94 


17.8x10.9 


@ 1 



Note. — Column 1 reports the maser species; column 2 gives the label number of the feature, as given in Tables [3] 
[4] and [5] column 3 lists the observing date; columns 4 and 5 report the absolute position; column 6 gives the peak 
brightness; column 7 reports the (major and minor) HPBW size and PA of the naturally-weighted (NW) beam. The 
PA of the beam is defined East of North. 
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Table 3. Parameters of VLB A 22.2 GHz H 2 maser features in G28.87+0.07. 



Feature epochs of Vlsr S„ Act AS Y a ~V$ 

Number detection (km s -1 ) (Jy) (mas) (mas) (km s _1 ) (km s _1 ) 






1,2,3,4 


103.17 




66.83±0.05 


0.005±0.05 








1 


1,2,3,4 


106.01 


53.99 








-16.9±1.0 


1.4±1.0 


2 


1,2,3,4 


105.07 


26.07 


37.44±0.05 


-3.02±0.05 


-23.0±1.2 


-26.4±1.2 


3 


1 


104.36 


25.38 


-2.88±0.05 


2.60±0.05 






4 


1,2,3 


105.27 


19.92 


-0.85±0.05 


0.25±0.05 


-2.8±2.2 


0.2±2.3 


5 


2,3,4 


104.27 


17.62 


-3.29±0.05 


2.79±0.05 


-29.8±1.7 


0.2±1.8 



Note. — Column 1 gives the feature label number; column 2 lists the observing epochs at which the 
feature was detected; columns 3 and 4 report the value of the intensity-weighted LSR velocity and flux 
density of the strongest spot, averaged over the observing epochs; columns 5 and 6 give the position 
offsets (with the associated errors) along the R.A. and Dec. axes relative to the feature #1, measured at 
the first epoch of detection; columns 7 and 8 give the components of the relative proper motion (with the 
associated errors) along the R.A. and Dec. axes, measured with respect to the reference feature #0. 
This table is published in its entirety in the electronic edition of the Astrophysical Journal. A portion is 
shown here for guidance regarding its form and content. 



Table 4. Parameters of EVN 6.7 GHz CH 3 OH maser features in G28.87+0.07. 



Feature 


epochs of 


Vlsr 


Fpcak 


I 


Aa 


AS 




vs 




detection 


(km s- 1 ) 


(Jy beam -1 ) 


(Jy) 


(mas) 


(mas) 


(km s- 1 ) 


(km s- 1 ) 


1 


1,3,4 


104.88 


0.155 


0.190 








1.7±6.7 


-15.6±8.1 


2 


3,4 


105.21 


0.268 


0.283 


-12.41±0.18 


-11.24±0.18 







Note. — Column 1 gives the feature label number; column 2 lists the observing epochs at which the feature was detected; 
columns 3 and 4 provide the intensity- weighted LSR velocity and the flux density of the strongest spot, respectively, averaged 
over the observing epochs; columns 5 and 6 give the position offsets (with the associated errors) along the R.A. and Dec. axes, 
relative to the feature #1; columns 7 and 8 give the components of the relative proper motion (with the associated errors) along 
the R.A. and Dec. axes, measured with respect to the reference feature #0 (the "center of motion") of water masers. 
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Table 5. Parameters of VLBA 1.665 GHz OH maser features in G28.87+0.07. 



Feature Pol V LSR F pcak I Ax Ay 







(km s x ) 


(Jy beam x ) 


(Jy) 


(mas) 


(mas) 


1 


L 


102.9 


1.94 


5.26 








2 


L 


106.1 


1.44 


3.52 


-233.57±0.11 


143.92±0.13 


3 


L 


105.4 


0.67 


1.78 


-139.46±0.17 


46.91±0.21 



Note. — Column 1 gives the feature label number; column 2 indicates the circular 
polarization of the maser emission; columns 3, 4, and 5 report the intensity-weighted 
LSR velocity, and the intensity and the integral flux of the strongest spot, respec- 
tively; columns 6 and 7 list the position offsets (with the associated errors) along 
the R.A. and Dec. axes, relative to the feature #1. 
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Table 6. Flux densities measured towards the G28.87+0.07 region. 



A 


F v 


instrument 


(cm) 


(Jy) 




3.6 


0.9 


Spitzer/GLIMPSE 


4.5 


2.1 


Spitzer/GLIMPSE 


5.8 


5.3 


Spitzer/GLIMPSE 


8.0 


9.4 


Spitzer/GLIMPSE 


8.3 


3.3 


MSX 


12 


8.6 a 


IRAS 


12.1 


7.1 


MSX 


14.7 


15.3 


MSX 


21.3 


34.0 


MSX 


24 


>21 


Spitzer/MIPSGAL 


24.6 


87.4 


Subaru/COMICS 


25 


105 a 


IRAS 


60 


1720 a 


IRAS 


70 


1790 


Herschel/Hi-GAL 


100 


<3180 


IRAS 


160 


1350 


Herschel/Hi-GAL 


250 


630 


Herschel/Hi-GAL 


350 


200 


Herschel/Hi-GAL 


500 


53 


Herschel/Hi-GAL 


870 


8.8 


APEX/ATLASGAL 


1100 


2.5 


CSO/BGPS 



a all IRAS fluxes are assumed to be upper limits in Figure [8] 
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7. ONLINE-ONLY MATERIALS 
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Table 7. Parameters of VLB A 22.2 GHz H 2 maser features in G28.87+0.07. (Table E] in 

the print version.) 



Feature epochs of Vlsr 




Aa 


A<5 


V Q 


v 5 


Number detection (km s -1 ) 


(Jy) 


(mas) 


(mas) 


(km s" 1 ) 


(km s _1 ) 






1,2,3,4 


103.17 




66.83±0.05 


0.005±0.05 








1 


1,2,3,4 


106.01 


53.99 








-16.9±1.0 


1.4±1.0 


2 


1,2,3,4 


105.07 


26.07 


37.44±0.05 


-3.02±0.05 


-23. Oil. 2 


-26.4±1.2 


3 


1 


104.36 


25.38 


-2.88±0.05 


2.60±0.05 






i 


1,2,3 


105.27 


19.92 


-0.85±0.05 


0.25±0.05 


-2.8±2.2 


0.2±2.3 


5 


2,3,4 


104.27 


17.62 


-3.29±0.05 


2.79±0.05 


-29.8±1.7 


0.2±1.8 


6 


2 


105.06 


15.08 


2.36±0.05 


-3.68±0.05 






7 


2 


105.14 


14.36 


14.86±0.05 


-146.72±0.05 






8 


1,2,3,4 


105.21 


12.16 


36.29±0.05 


-3.33±0.05 


0.5±1.3 


-10.0±1.3 


9 


1,2,3,4 


105.45 


11.21 


-1.71±0.05 


0.79±0.05 


-26.3±1.5 


7.0±1.5 


10 


1,2,3,4 


105.25 


9.40 


21.35±0.05 


-2.95±0.05 


0.9±1.5 


-5.7±1.5 


11 


1 


103.92 


9.39 


-3.79±0.05 


3.17±0.05 






12 


1,2,3,4 


104.80 


9.30 


29.94±0.05 


-5.50±0.05 


8.2±1.5 


-6.9±1.6 


13 


1,2,3 


104.00 


9.15 


-5.85±0.05 


4.22±0.05 


-21.5±2.2 


-13.5±2.3 


14 


2 


105.35 


8.57 


-1.00±0.05 


24.53±0.05 






15 


1,2,3,4 


106.34 


8.22 


56.56±0.05 


-3.51±0.05 


-21.2±1.4 


3. Oil. 5 


16 


1,2,3,4 


106.73 


7.53 


-1.84±0.05 


0.91±0.05 


-25.3±1.1 


6.4±1.2 


17 


2 


105.12 


7.33 


3.31 ±0.05 


-46.40±0.05 






18 


2,3,4 


104.94 


7.23 


29.24±0.05 


~165.12±0.05 






19 


1,2,3,4 


105.87 


7.22 


34.93±0.05 


-5.68±0.05 


21.1±1.6 


-25.0±1.6 


20 


2 


104.70 


7.02 


-11.90±0.05 


192.10±0.05 






21 


1,2,3,4 


105.31 


6.90 


29.01±0.05 


-5.21±0.06 


7.6±1.5 


-0.7±1.6 


22 


1,2 


104.23 


6.82 


36.56±0.05 


-3.49±0.05 






23 


1,2,3,4 


105.42 


6.57 


27.79±0.05 


~4.98±0.06 


8.3±1.6 


-12.6±1.7 


21 


2 


105.11 


6.48 


-0.06±0.05 


-42.00±0.05 






25 


2 


105.10 


6.45 


15.54±0.05 


-104.30±0.05 






2(5 


1,2,3 


95.16 


6.20 


202.71±0.05 


50.48±0.05 


55.3±1.8 


18.4±1.9 


27 


1,2,3,4 


105.50 


6.16 


21.32±0.05 


21.04±0.05 


11.5±1.9 


-11.6±1.9 


28 


2 


104.36 


5.79 


-1.96±0.05 


2.10±0.05 






29 


1,2,3,4 


103.91 


5.57 


-6.62±0.05 


4.52±0.05 


-42.5±1.3 


-4.2±1.3 



Table 7— Continued 



Feature epochs of Vlsr Si/ Act AS ~V a ~Vs 

Number detection (km s _1 ) (Jy) (mas) (mas) (km s _1 ) (km s _1 ) 



30 


2,3,4 


103.78 


5.30 


-2.48±0.05 


2.28±0.05 


4.2±1.9 


-23.6±2.0 


31 


1,2,3,4 


105.13 


5.28 


31.91±0.05 


-9.86±0.05 


-5.2±1.8 


-26.4±2.0 


32 


1,2,4 


104.22 


5.23 


-2.89±0.05 


-40.09±0.05 


-50.6±2.7 


-8.5±2.7 


33 


1 


104.35 


5.23 


-8.60±0.05 


8.75±0.05 






34 


1,2,3,4 


106.26 


5.21 


21.33±0.05 


-4.99±0.05 


-1.6±1.6 


-6.3±1.6 


35 


3 


104.71 


5.17 


14.22±0.05 


-239.02±0.05 






36 


1,2,3,4 


105.46 


4.86 


20.48±0.05 


-4.65±0.05 


4.1±1.7 


-3.7±1.7 


37 


1,2,3 


105.15 


4.78 


28.10±0.05 


11.48±0.05 


1.9±3.1 


-8.7±3.2 


38 


2 


105.54 


4.72 


36.95±0.05 


13.71±0.05 






39 


3 


104.72 


4.57 


-10.01±0.05 


121.21±0.05 






40 


1,2,3,4 


108.62 


4.46 


2.52±0.05 


-2.55±0.05 


-7.9±1.2 


-6.2±1.2 


41 


2 


104.93 


4.37 


36.76±0.05 


-27.47±0.05 






42 


1,2 


104.71 


4.32 


-5.05±0.05 


46.44±0.05 






43 


2 


105.12 


4.17 


25.61±0.05 


0.01±0.05 






14 


1 


104.23 


4.06 


-6.09±0.05 


30.67±0.06 






45 


1,2,3,4 


106.20 


4.03 


61.07±0.05 


-7.08±0.05 


6.2±1.7 


-5.8±1.7 


46 


2 


105.52 


3.96 


33.94±0.05 


43.85±0.05 






47 


1,2,3,4 


104.67 


3.95 


-2.36±0.05 


26.39±0.05 


-50.8±1.8 


13.5±1.8 


48 


1,2,3 


104.04 


3.91 


-4.62±0.05 


3.58±0.05 


-41.3±2.5 


7.8±2.7 


49 


1,2,3,4 


103.85 


3.84 


-6.13±0.05 


47.09±0.05 


-31.9±1.7 


10.5±1.7 


50 


1,2,3,4 


105.44 


3.71 


21.97±0.05 


-20.76±0.05 


1.2±1.9 


1.5±2.1 


51 


1,2 


105.06 


3.67 


38.10±0.05 


-20.04±0.05 






52 


1 


105.41 


3.60 


-0.82±0.05 


-23.65±0.05 






53 


2,3,4 


105.55 


3.24 


35.42±0.05 


-29.60±0.05 


10.5±2.7 


-22.8±2.7 


54 


1,2,3,4 


106.82 


3.22 


23.89±0.05 


-3.16±0.05 


-3.5±1.4 


-9.8±1.4 


55 


2,3 


106.46 


3.13 


13.73±0.05 


-123.75±0.05 






56 


4 


105.35 


3.07 


3.11±0.05 


-47.13±0.06 






57 


3,4 


105.95 


3.02 


-4.34±0.05 


5.38±0.06 






58 


2,3,4 


106.37 


2.93 


4.53±0.05 


-2.16±0.05 


-38.8±2.4 


-21.0±2.4 


59 


3 


105.39 


2.90 


25.72±0.05 


-25.59±0.05 







Table 7— Continued 



Feature epochs of Vlsr Si/ Act A<5 V a Vs 

Number detection (km s _1 ) (Jy) (mas) (mas) (km s _1 ) (km s _1 ) 



60 


2 


105.01 


2.86 


0.23±0.05 


1.02±0.05 






61 


1 


104.90 


2.81 


-1.77±0.05 


2.26±0.05 






62 


4 


105.34 


2.77 


2.65±0.05 


-5.13±0.06 






63 


1,2,3,4 


105.46 


2.72 


37.26±0.05 


20.82±0.05 


-9.8±1.8 


-12.8±1.S 


64 


3 


105.83 


2.63 


35.05±0.05 


-48.30±0.06 






65 


1,2 


105.30 


2.51 


57.47±0.05 


-3.62±0.05 






66 


4 


103.22 


2.40 


-5.42±0.05 


3.83±0.05 






67 


1,2,3 


105.21 


2.34 


38.23±0.05 


-3.01±0.05 


-14.0±3.5 


12.2±3.7 


68 


1,2,3 


106.13 


2.30 


20.32±0.05 


-47.78±0.05 


-10.4±3.1 


2.4±3.3 


69 


1,2,3,4 


106.85 


2.28 


3.34±0.05 


-2.71±0.05 


-3.2±1.4 


-3.8±1.5 


70 


2,3 


105.04 


2.26 


40.16±0.05 


-5.89±0.05 






71 


2 


106.54 


2.21 


26.23±0.05 


-8.07±0.05 






72 


4 


105.71 


2.15 


0.42±0.06 


-42.97±0.06 






73 


3 


105.30 


2.08 


6.51±0.05 


-1.41±0.05 






74 


1,2,3,4 


106.74 


1.97 


22.28±0.05 


-2.91±0.05 


-5.8±1.7 


0.8±1.7 


75 


2 


106.40 


1.95 


56.00±0.05 


-3.39±0.05 






76 


3 


106.37 


1.91 


20.16±0.05 


-23.71±0.06 






77 


1,2,3,4 


104.43 


1.87 


23.93±0.05 


-4.51±0.05 


-9.5±1.8 


-2.7±2.0 


78 


3 


106.05 


1.87 


-4.11±0.05 


71.74±0.06 






79 


3 


106.54 


1.82 


8.25±0.05 


35.04±0.05 






80 


3 


106.27 


1.82 


20.56±0.05 


-65.42±0.05 






81 


3 


105.94 


1.81 


-1.45±0.05 


43.02±0.06 






82 


2,3 


106.59 


1.80 


3.66±0.05 


-45.51±0.05 






83 


1,2,3,4 


94.90 


1.80 


138.95±0.05 


67.32±0.05 


6.7±1.2 


14.6±1.2 


84 


3 


106.57 


1.72 


23.19±0.05 


-46.20±0.05 






85 


4 


105.12 


1.67 


38.38±0.06 


33.81±0.06 






86 


1 


105.73 


1.66 


22.67±0.05 


-4.96±0.05 






87 


3 


106.06 


1.66 


2.01±0.05 


-1.62±0.06 






88 


3 


106.37 


1.64 


20.18±0.05 


1.15±0.05 






89 


4 


105.42 


1.61 


29.71±0.06 


-98.11±0.06 







Table 7— Continued 



Feature epochs of Vlsr Si/ Act AS ~V a ~Vs 

Number detection (km s _1 ) (Jy) (mas) (mas) (km s _1 ) (km s _1 ) 



90 


1,2,3,4 


95.86 


1.60 


220.77±0.05 


55.75±0.06 


53.6±1.4 


17.8±1.5 


91 


1 


103.89 


1.58 


0.44±0.05 


-2.15±0.05 






92 


1,2 


104.29 


1.52 


23.52±0.05 


-4.52±0.05 






93 


4 


105.12 


1.49 


31.87±0.06 


-34.73±0.06 






94 


1,2,3 


104.78 


1.48 


7.67±0.05 


-1.02±0.05 


-13.0±3.9 


-12.2±4.1 


95 


2,3,4 


107.30 


1.45 


-3.86±0.05 


1.76±0.06 


-27.6±2.3 


7.0±2.4 


96 


4 


105.31 


1.39 


29.71±0.06 


-55.31±0.06 






97 


1,2,3,4 


106.59 


1.37 


5.76±0.05 


-2.52±0.05 


-22.6±1.8 


-9.8±1.9 


98 


3 


103.95 


1.35 


-6.91±0.05 


-38.75±0.06 






99 


2,3 


106.78 


1.30 


46.15±0.05 


0.89±0.05 






100 


4 


104.94 


1.28 


30.01±0.06 


-31.50±0.06 






101 


2 


106.96 


1.25 


52.39±0.05 


-8.39±0.05 






102 


1 


103.88 


1.25 


-6.46±0.05 


7.37±0.07 






103 


2,3 


106.11 


1.24 


-0.24±0.05 


3.29±0.06 






104 


3 


104.07 


1.12 


-0.25±0.05 


-44.67±0.06 






105 


4 


105.11 


1.11 


33.36±0.06 


2.56±0.06 






106 


4 


105.50 


1.11 


37.49±0.06 


-29.30±0.06 






107 


2,3 


111.99 


1.10 


71.41±0.05 


0.79±0.06 






108 


1,2,3,4 


106.79 


1.09 


49.43±0.05 


-3.39±0.05 


6.3±1.7 


-14.1±1.7 


109 


2,3 


108.68 


1.03 


100.94±0.05 


20.38±0.06 






110 


3 


106.57 


0.97 


1.57±0.05 


-46.93±0.05 






111 


1 


106.80 


0.94 


1.44±0.05 


-1.20±0.05 






112 


1,2,3 


106.82 


0.86 


48.34±0.05 


-3.64±0.05 


2.5±4.0 


-9.2±4.2 


113 


4 


106.91 


0.83 


-1.46±0.06 


-42.03±0.06 






114 


3 


104.67 


0.81 


23.33±0.05 


-4.89±0.06 






115 


4 


103.55 


0.79 


-11.43±0.05 


10.12±0.06 






116 


4 


103.92 


0.73 


-7.58±0.06 


8.39±0.06 






117 


1,2,3,4 


103.66 


0.72 


179.75±0.05 


-103.20±0.05 


18.9±2.2 


1.3±2.4 


118 


1,2 


94.32 


0.69 


139.93±0.05 


66.86±0.05 






119 


4 


103.64 


0.68 


0.34±0.06 


-114.27±0.06 







Table 7— Continued 



Feature epochs of Vlsr Si/ Act AS ~V a ~V& 

Number detection (km s _1 ) (Jy) (mas) (mas) (km s _1 ) (km s _1 ) 



120 


3,4 


107.47 


0.66 


167.34±0.05 


-261.03±0.06 






121 


1,2,3 


107.35 


0.61 


50.20±0.05 


-3.27±0.05 


-0.6±3.7 


-9.6±4.0 


122 


1 


107.18 


0.57 


-3.75±0.05 


0.95±0.05 






123 


1,2,3 


94.48 


0.57 


214.01±0.05 


54.68±0.06 


48.1±2.6 


18.7±2.8 


124 


3,4 


107.11 


0.53 


-6.35±0.05 


6.41±0.07 






125 


2 


95.24 


0.49 


202.75±0.05 


51.60±0.05 






126 


1 


107.01 


0.48 


-5.96±0.05 


6.03±0.05 






127 


1,2,3,4 


94.75 


0.41 


206.85±0.05 


51.60±0.06 


43.1±1.5 


14.0±1.6 


128 


4 


108.82 


0.38 


-1.59±0.06 


3.18±0.06 






129 


1,2,3,4 


93.64 


0.38 


140.63±0.05 


66.57±0.05 


5.9±1.6 


14.3±1.6 


130 


3 


103.41 


0.38 


-6.86±0.06 


74.45±0.07 






131 


2,3,4 


107.11 


0.37 


84.03±0.05 


-0.53±0.06 


-0.8±2.7 


-15.5±2.9 


132 


2 


95.94 


0.36 


138.45±0.05 


67.75±0.06 






133 


1,2,3 


93.60 


0.36 


213.24±0.05 


55.03±0.11 


46.1±2.9 


21.3±3.1 


134 


1,2 


99.32 


0.34 


121.98±0.05 


91.75±0.06 






135 


1,2,3,4 


93.62 


0.34 


205.00±0.05 


51.40±0.06 


42. Oil. 7 


11.4±1.7 


136 


1,2,3,4 


95.47 


0.34 


204.63±0.06 


51.38±0.06 


52.1±1.8 


17.2±2.1 


137 


1,2,3 


109.59 


0.33 


98.29±0.05 


11.15±0.06 


-12.3±3.0 


-12.6±3.4 


138 


3,4 


99.16 


0.31 


122.29±0.05 


91.99±0.06 






139 


4 


93.97 


0.31 


214.80±0.06 


55.26±0.06 






140 


3 


103.39 


0.30 


58.58±0.06 


-3.38±0.07 






141 


3 


107.37 


0.28 


65.53±0.05 


5.02±0.06 






142 


1,2 


95.14 


0.27 


202.97±0.05 


49.48±0.05 






143 


1,2,3 


107.23 


0.25 


166.00±0.05 


-257.48±0.06 






144 


1,2,3 


108.39 


0.25 


45.92±0.05 


-6.11±0.05 


2.7±3.7 


-8.6±4.4 


145 


4 


95.90 


0.25 


218.01±0.06 


61.98±0.07 






146 


1,2,3 


94.16 


0.23 


136.39±0.05 


74.76±0.06 


6.9±1.9 


7.9±2.2 


147 


1,2,3 


111.87 


0.22 


171.96±0.05 


-42.98±0.05 


20.3±2.3 


-6.5±2.6 


148 


1,2 


99.12 


0.22 


122.43±0.05 


91.64±0.05 






149 


4 


107.41 


0.21 


-8.25±0.06 


7.57±0.07 







Table 7 — Continued 



Feature epochs of Vlsr S„ Ack A<5 V q Vg 

Number detection (km s _1 ) (Jy) (mas) (mas) (km s _1 ) (km s _1 ) 



1 *,n 

10U 


A 
I 


lUo.zo 


u. 


01 
Z 1 


a no-uo H7 
— O.UUitU.U / 


oo oc_i_n n7 
— ZZ. ZoztU.U/ 


1 -1 1 
1 u 1 


9 


yo.zu 


u. 


"1 


i aq-uo o^ 


11 1 /I— 1-0 07 


152 


4 


108.66 


0. 


19 


2.04±0.07 


-4.78±0.07 


153 


1.2 


102.57 


0. 


16 


100.89±0.05 


8.93±0.05 


154 


1,2,3 


95.22 


0. 


16 


202.79±0.06 


56.68±0.05 


155 


2,3,4 


93.33 


0. 


16 


205.03±0.06 


58.09±0.06 


156 


1 


99.90 


0. 


14 


99.55±0.05 


9.10±0.06 


157 


4 


103.02 


0. 


14 


-2.92±0.07 


-0.52±0.08 


158 


1,2,3 


96.62 


0. 


13 


217.56±0.05 


54.15±0.06 


159 


1 


95.29 


0. 


13 


202.24±0.06 


47.80±0.06 


160 


2 


107.64 


0. 


13 


182.83±0.05 


-174.07±0.07 


161 


3 


108.71 


0. 


13 


2.47±0.06 


-4.01±0.05 


162 


2 


107.42 


0. 


12 


170.05±0.06 


-263.08±0.08 


163 


4 


94.87 


0. 


12 


135.39±0.06 


73.40±0.08 


164 


2 


94.98 


0. 


12 


201.51±0.06 


43.74±0.08 


165 


3 


108.47 


0. 


12 


82.90±0.06 


1.27±0.09 


166 


3 


95.84 


0. 


11 


225.25±0.06 


7.28±0.10 


167 


2 


98.34 


0. 


11 


128.73±0.06 


86.69±0.09 


168 


3 


107.50 


0. 


10 


-1.99±0.07 


2.06±0.06 



46.4±3.0 
27.8±4.5 



51.8±2.0 



1.8±3.7 
-1.8±4.5 



18.0±2.4 



Note. — Table |3] in printer. Column 1 gives the feature label number; column 2 lists the observing 
epochs at which the feature was detected; columns 3 and 4 report the value of the intensity-weighted LSR 
velocity and flux density of the strongest spot, averaged over the observing epochs; columns 5 and 6 give 
the position offsets (with the associated errors) along the R.A. and Dec. axes relative to the feature #1, 
measured at the first epoch of detection; columns 7 and 8 give the components of the relative proper motion 
(with the associated errors) along the R.A. and Dec. axes, measured with respect to the reference feature 
#0. 



